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On the Variation of Uncanonical Arbitrary Constants; with an 
application to the Planetary Theory. By A. Y. G. Campbell, 
B.A., Scholar of Trinity College, Cambridge. 

(Communicated hy Sir li. S. Ball.) 


In attacking the planetary or lunar theories by the method 
of the variation of arbitrary constants, it is necessary to evaluate 
several functions usually denoted by [a, e], &c., unless a canonical 
set of arbitrary constants is chosen. In this paper I endeavour 
to show that the equations giving the variations of these con¬ 
stants can be found more quickly without evaluating the 
functions \a, e], &c. It is necessary to find one function 
(denoted in this paper by C with a suffix) corresponding to 
each arbitrary constant. In the first two sections the theory 
and its application to the elliptic constants of the planetary 
theory are briefly explained. In the following section it is 
shown that if any solution of the new equations is found we can 
proceed to find the equations, giving the variation of the 
new arbitrary constants. The next approximation can be ob¬ 
tained in a manner similar to that used in the previous case. 
It is shown also in the last section that in certain cases there 
exists a method of solution similar to Jacobi’s method of solution 
of the Hamiltonian canonical equations. Among these eases is 
included the set of equations giving the variations of the elliptic 
elements in the planetary theory. 


i. 


Let 


«1J a 2) 


solution of the canonical equations 


a 2 n be the set of arbitrary constants of any 


. dH . dH,, 

9i dp.' Pi dq 1 ‘ 2 ‘ 


0) 


It is required to find the equations giving the i ates of Variation 
Of a,, a 2 , . . . a 2re when a disturbing function O is added to H. 

If L is the Lagrangian function corresponding to H, 

dL . dL ■ , ■ . 

and if a r is one of the arbitrary constants in terms of which L is 
supposed to be expressed, 

dL f-f /dL dq i ^ dL dq\\ 
da. i=1 \dq i da r dq. da) 


r % 
i=n 


dt , 1 da / 
1=2:1 
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Jan. 1897. Uncanonical Arbitrary Constants. 119 

the terms in the bracket being expressed in terms of the arbitrary 
constants, &e. In the above q { have been written for the 
values of 

d 'Pj, *l t> 

dt dt 

when p i} q t are expressed in terms of a 1? a 2? . • . ci 2m t. 

It follows that 


IteSK 


dL 

da 


dt + A r 


dL 

cla r 


being expressed in terms of a u a 2 , . . . a 2n , t, and integrated 
with regard to t explicitly ; A r is a function of a,, a 2 , . . . a 2n) and 


a l> ®2> ‘ ‘ • 

but must not contain t. 

Now 

® a r J ti i=ti 

if we put 

A r +(L) fr+B, = o 

t=t x da r 

we have 


and similarly, 

X(pA) = d - f lik-b, 

l daJ da J t. 


dty 

da,. 


( 2 ) 


If we differentiate the first of these expressions with respect to 
ct s and the second with respect to a r) on the supposition that all 
the functions involved are expressed in terms of a l5 a 2 , . . . a 2m , t, 
we find that 



dB dL 

S _ r 

da r da g 


( 3 ) 


where 


[®r» ®»] : 


X( dqi dp i) 

. 'da r da s da s da r J 


( 4 ) 


Since B r , B s are independent of the time explicitly, it is clear 
that [a„ a s ] does not contain t explicitly when it is expressed in 
terms of a 1} a 2 , . . . a 2n . It follows also that we can assign any 
value we please to t in equation (2) as far as regards the deter¬ 
mination of B r 
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If 

.«> 

where A is any arbitrary function of a l5 a 2 , . . . a 2 „ it is clear 
that 


= dCs dC r 
da r da s 


( 6 ) 


If we put t=t l in equation (2) 



• ( 7 ) 


Having thus determined [a r , a s ], we can now proceed in the 
usual manner to find the new equations. 

Let the original equations be 


dpi dqi 


. . ( 8 ) 


Then, since a,, a 2 , . . . a Un have been so chosen that when 
q ly q 2 , ■ . . q n , P\, p2) • • • Pn are expressed in terms of them and t, 


dxp,_dH. dpi 

dt dpi dt 


dK 

dqi 


(*= I, 2, . . . n) 


it follows that 


da r ' dpi da r U ' dos 

r=l r=l 1 

Multiplying each of these equations by 


(i=i, 2, 


. . n). 


dpi dqj, 
dal da 


respectively, and adding the complete set, we see that 


dSl 

da s 


■r—2n 

= 2 Or, 


«s]«r 



are the required equations. 
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§2. Application to the Planetary Theory. 

The motion of the planets is governed in general principally 
by the. attractive force of the Sun, which would, if acting alone, 
cause each planet to describe an ellipse. The deviations from 
the elliptic orbit are small; it is an advantage to deal with 
variables whose rates of variation are small, for it is then possible 
to use methods of approximation. In some expositions of the 
planetary theory the first change of variables is to the constants 
of the elliptic orbit. The set which are usually denoted by a , e, e, 
g, a, <f>, are uncanonical. 

Let the origin of coordinates be taken at the Sun’s centre, 
which is supposed to be at rest or moving with uniform 
velocity in a straight line. If (x, y, 2) are the coordinates of the 
planet referred to three fixed rectangular axes through the 
origin 

L = 7 (i- + y- 4 -i 3 ) ~- 

2 r 

when the units are properly chosen. 

Then 

Pl = X,p 2 = y,p 3 = s 

H =*(?> ! 2 +Jv+P 3 2 ) 

2 r 


Let b denote any one of the constants a , e, e, g, a, <£. It is neces¬ 
sary to find 




dz 


db db 

i 

We shall now consider the transformation to new rectangular 
axes— 

x = + l 2 v + hC | 

y = m 1 l- + m 2 7i + m 3 ( h 

z = + n 2 t) 4 - n s C / 

where the direction cosines may be functions of the elements, 
but not of t explicitly. 


»dx , • lit! 
J - 1 - 


db'* db 
+ 


• / 7 > 

±* s .(ht- 


+ 

:d£ t .dr) 


+ i« + u)(*,g + '.g + 4 f + 4 + 4 + f|) 


dC 


*db + V db + ^db 


r dl x dm. , „„ dn,\ 
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Let the plane x y cut the plane of the ecliptic in the axis £, 
and let the axis 17 be in the plane of the ecliptic. The angle 
between x and | 



is a, and if X denote the axis which is taken along the major 
axis of the ellipse, the angle between £ and X is g ; the inclina¬ 
tion of the ecliptic to the plane x y is cj>; a and e have their 
usual meanings in an ellipse, and e is the epoch. 

The solution of the elliptic motion is given by the following 
equations— 

X = «(cosw — e) \ 

Y = #a/i —e 2 . sin u r 
a~%t+e = u — e sin u ' 

| = X cosy—Y siny ) 

V = X sin g + Y cos g \ 
x — % cos a — 7j sin a cos <p * 
y = £ sin a + k\ cos a cos (f> - 
# = 7] sin (p j 


Considering the description of area we see that 

^ = XY—YX= 


Hence 


. (lx . (1)! . ds 

dtp dtp dtp 

■ dx.dy . dz _ 

=f §+«t-XY-YX 
off dg dg d£ dg 


= a /ct (1 — e 2 ) 


and when b—a , e or e 


+ + -ir dX -,y- dY 

* db + y db + *db~*~db + Y W 
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We put 
Then 


Uncanoniccd Arbitrary Constants. 


L = H+ z = 2 —L. 
r r 2 a 


t l = — eat, so that tt-o when f ==t l 




a(i — - 2 a 


Hence using the formula 


rfA. 


?*=1 1 

we see that, if we assume A to be independent of g , a, <f>, 


C g — — Va (1 — e 2 ); G a = — Ja (i — e-) . cos ; 0 ^ = O. 

We have still to find C a , C c , C c . 

Since 

X=a (cos u—e), 


Y = aVi — & sin u, 
a~U+ e = u—e sin u 
(X) M=0 = a( i-e), (Y) k=0 = o 

/<*X\ _ n tdY\ 

\du ) u=0 ’ \du) u =°~ a '' /l ~ 

(X)»= 0 =■ o, and 




Hence 




„ / 2 1 \ dt. _ /^zA 

V*(i-e) 2 a/ + O \d6 / ), 1=3 + 

«=*■ 

C a =3 € _L/_ *__i\ 3 I .±i’ + ^4 = 3_! = + dA - 

* 2 ^\l— e 2/ 2 V a 1— e da 4 <£a" 


C,= 


«!A 


_/ 2 l\ 1 + /’ , rfA v 7 ® dA 

Ce = -v/fl ( - ) — a/ <2- t" “T- H -T" • 

\l — e 2 / 1 — p tfe 2 <£e 

If we put A= — \ '/». €, we see that 

C a =A-A C« = o, Ce = o, 

2 Va 


C 9 = — ■/«(! — e 2 ), C a = — Va(i — e 2 ) cos <f>, C^ = o. 


125 
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Hence from equation (9), the new equations are 

dn 1 . 1 , _ . 

— = -a 2e + 2 a * v'i —e-(^ + d cos <f>) 


da _ 

de 

da 1 1. 

— zz-a~m 
de 2 


- ah- 


a/I* 


: (g + a cos $) 


da I ,-- e ; 

— = —ar*V 1 — e 2 .a + a* • e 

dg 2 Vi—e z 

da 1 , --* . 6 

~T = — a~*V 1 —e 2 . a cos <p + a$—=— 
da 2 J\—e 2 


LYIL 3, 


e cos <f> 4- Va(i—e 2 ) sin <£ 


d(p 


= — a/& (1 — e 2 ) . sin $ . a. 


§ 3. Method of Successive Approximation . 

In proceeding to consider consecutive terms of the disturbing 
function, it is assumed that a solution containing 2n arbitrary 
constants /3 ]5 £ 2 , . . fi 2n has been found of the equations 


f^=C s - 2 ^,( 5 = I, 2, ..2«> 

r=l 


(10) 


and it is required to find the equations for y 3 b /? 2 , . . . /? 2 „ when 
0 2 is added to O,. 

Let 

r=2n. 

— Xij -4~ C r ct r . • . . • • . . (11) 

r=l 

dL, da s r ^°^ n dCr . das + r ^ c 

40 » * cU s dH k * ^ da, r dp k * r dp k 

S ~1 f=l S=1 »’*=! 



the terms in the bracket being expressed in terms of j8 b /S 2 , . .. 
£ 2lt , t, and C„ a, being written in the above for the values of 


dC s dag 
dt J dt 


when C w a s are expressed in teyms of / 3 1? /? 2 j • • m ^ 
It follows that 


s=I 
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I2 5 


—1 being expressed in terms of / 3 1? fi 2 , ■ ■ ■ ft 2n , t, and integrated 

dji k 

with regard to t explicitly ; D fc is a function of / 3 lf /? 2 , . .. ( 3 2n , 
but does not contain t explicitly ; t x may be any function of 
. . . / 3 2 »> hut must not contain t. 

If 

IV + (L,) ( _. — ' - E. 


d& : . 


then 

similarly 

Hence 

where 




5=1 

s~2n 


A_ d j 

■t \ 

[ L ,df — E/, 

! f. 

J rf/BtJ 

V 

■ \ d \ 

[ * L, dt — K / 

If) dfy 

1 * ~ 

1 

/ 


da , 


r a Q-\ _dEj 

&] (/3 (/3 


( 12 ) 


“ W; c/j 3 * 

5=1 


(13) 


(14) 


Since E ; , E ft are independent of t explicitly, it follows that 
[/3,-, /?,] does not contain t explicitly, when expressed in terms of 
/? 1} y8 2 , . . . y8 2w . It follows, also, that, as far as regards the de¬ 
termination of E; £ , we can assign any value we please to t in 
equation (12). 

If 

f *" e ‘ + £.<«> 


when B is any arbitrary function of fi u [i. 2 , . . . ft 2n , it is clear 
that 


[£V, ft] = 


dYj 


di\ 

dft 


and if we put t=t x in equation (12) 


*»=- 



dtf I 


(i6> 


(* 7 > 


Having thus determined [/?,,, /?,] we can now proceed to find the 
new equations. 

Let the original equations be 


d (Gj 4~ ^ 2 ) 
da$ 


=2 n 


. ^ a c r . 7 

— C s — ^ Of, (s — 1, . . 2 n) 


r —1 


da. 


. (18) 


L 
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Then, since ft, ft, . . . ft„ have been so chosen that when a 1( a 2 , 
. . . a> n , Cj, C 2 , . . . C 2 », are expressed in terms of them and t 

da , dC s ^ 

7«7 = dt da s ' W’ ^ = I> 2 ’ ’ ' ' 2n ) 

r —1 


it follows that 

dn 2 ^POC'-T R ^ r ~ 2n dCr da, . 
da s - f =i dfi^- * * da s ‘ dh • 


Multiplying this equation by 


das 

d&j 


(S=I, 2, 


. . 2 ») 


and adding the complete set we see that 


da. k ^ n s ^ n dC s da s 6 


Jc=%n s= 2 n 


"5 * 


^ ,l r y re ^C r aft 

Z r Z *% 


K —2 n 

= 2 [^A> £/] 

*=1 




lt= 2 n ijy 
k=l J 


2n) 


• • (19) 


are the required equations. 

It is clear that that this theorem includes that of §1 as a par¬ 
ticular case. If it is possible to express ft, ft, . . . ft„in terms of 
F„ F a , . . . F 2n , we see, by multiplying equation (19) by 

d&j 

d^m, 

and adding, that 

da, a d$ k , 

-JV' = ~ S • 4 («=t2... 2W) .... (20) 

k=l 

There is a species of reciprocity between the set of variables 
(Fj, F 2 , . . . F 2w ) and the set (ft, ft, . . . ft w ). 
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From a consideration of the coefficient of F, in (20), it is clear 
that 

nc 1 tt 1_ d fi r d&r 

^ K ' ] - ~dV ' ,/F,* 


§ 4 - K n of the functions denoted by C, say C, i+ ,, C n+ .,,.. . C 2n , 
are functions of the n «s only, which do not correspond to them, 
viz. a h a 2 , . . . a n , it is possible to indicate a method of-solution of 
the equations : 

A 

"T~ = ^ 
da s 

which is similar to Jacobi’s method of solution of the canonical 
equations. 

Let fi x , fi 2 , • • • fin, 7 u 725 •• • 7 h> be 2 n variables connected 
with aj, a 2 , . . . a 2n by the equations 


r—2n rp 

M»-l, 


(9) 


dij/ 

dai 

dxp 

d 6 i 


r—%n 


C, CL.f 


r=n+1 


dC r 

da; 


7i 


i 

f (< = I- 2, . 

) 


. . n) 


( 21 ) 


where ^is a function of a,, a 2 , . . . a„, fi u fin ■ ■ ■ fin, t only. 
Let S denote an arbitrary variation of the variables. 
Then since 

SG r = X a lfisa i+ dC fSf 

. . da L dt 

t=i 4 


if C, do not containa M+1 , a m+2 ,. .. o. 2 „, where r=n+ 1, n+ 2, 



d&i 



dt 


2 n , 


when 


i=n i=2n 

= 2 (Q + yAPi) — 2 } OjS Cj + <pst 

1=1 Wi + 1 



i=?i + l 


dC 


(22) 


Let A denote another independent arbitrary variation. 

Then 

i—n 

^ [AC^a^ -f A+ C^ASa* 4 - y^A^] 

£=1 

— ^ (Aa^O^ + cs^A^C^) 4* A<£ . dt 4- (ft , Adt 

i—nA-1 

Similarly 

i=n 

dA$ — ^ [dCiA&i + djiA&i 4 - C^Aa* 4- y i $A^] 

4=1 

i=2n 

— ^ [SojAOi + o^SAC^] + d<(> . At 4- <p , dAt. 

i*=n +1 

L 2 
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But 

A 8 tp = 8 Aip, 8 Aai = A 5 c^, &e. 


Hence 



i=2n 

r=/i 



[So*ACi - AaiSCi] + V[ - 8 y r A$ r + A 7r 8^.] + A<p .St — S(j).At = 0 . (23) 


i =1 i-l 

Now a,, a 2 , . . . o-2n can be expressed in terms of /8 1} /? 2 , . . . 
j3„, yi, y 2 , • • • 7n, t, and vice versa. 

Let ( )' denote that the quantity enclosed in the bracket is 

expressed in terms of / 3 2 , . . . ft n , y u Jp y n, t and let all quan¬ 
tities except if/, which are not enclosed in ( )', be in future ex¬ 

pressed in terms of a i} a 2 , . *. a 2m , t ; but let if/ be in all cases 
expressed in terms of a 1? a 2 , . . . a n , / 3 t /? 2 , . . . /?„, t. 


Then 



Equation (23) can then be written 




s—1 r =1 


i=* 1 


+ fHj 


=1 

■r~n 


^mm ASr+ ^ y)+ m. 


aw 


r—1 

78^1=°. (24) 

1 ~z x a&i dt ) 

But the variations 

8a„ Sa 2 , . . . 8a,„,, Si, Aft,, A 0 2 , . . . A/ 3 n , Ay,, Ay.,, . . . Ay,,, At 


are all arbitrary and independent; each coefficient of these 
quantities in equation (24) must therefore be zero. 
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Hence we find : 

W _ dC_dy r = Q \ 

d(3 r ^ d\ 3 r do., dai j 

6=1 ; 

djGiY _d(aj dC, _ dp r = o ■" 

dy r * dy r da, da, 

d(Ci)' S ^ n d(a s y dc,_d<p 

dt * 


5=1 


dt da, da . 


- t kh- 5 <*(«»)' - d w = o 

<ZZ ^ dp, ' rff rf/0 r 

d-P, _ s ^ d(a s )' dC, d(<p ' _ 
Z d 7r ' dt ‘ dy T 


5 1 


d(<p)>_^ d(a s Y dC,_dt 
dt dt ' dt dt 

5=1 


Now the original equations are 


Hence 


da A r=2n d Q 
r=l 

fiWr * L dP r ■ 3 A eft ‘ ^7 • J 

s=l »=1 

*^"p*C#.y * -1 dtCiY 

5=1 - - 


i=l 


* _d(c s y A -| 

“^ L^87 ,0s ^sr- as J 

5=1 

= ^ r dC, d(a s y *=*» ^G Si (Z(a,)' 

2* \_dt ‘ Zf <Za £ s dp r 

s—1 %—1 

/ s === 2?& 1 / \ / 7 .n ^ - 

= -*b+m-'%A by( 2 7) 


dt dfir ** v da{ 

i =1 


whence we obtain 
Similarly 


>=-—(4»-ay 1 

1 dp r - ’ 

Pr^J-^-ay 
dy r 


129 


■ (25) 


(26) 


(27) 


(28) 


(9) 


(»•=»= 1, 2 , ... n). . (29) 


These new equations are canonical. 

If <£=0, we see that ft, ft, . . . ft, y u y 2 , . . . y K , form a 
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canonical set of arbitrary constants, the solution being given by 
the equations 



- (i= 1,2, ... n) . . . . (21) 


when 



(22) 


i—n + 1 


From the equations 



we can in general find a n+l , a n+2 , . . a 2n in terms of 



If these values are substituted in the equation 



(30) 


we obtain a differential equation of the first order, for if/ in terms 
of a,, a 2 , . . . a„, t. If any complete solution of this equation, con¬ 
taining n arbitrary constant /?,, /3 2 , . ■ . be found, the solution 
of the equations (9) is given by equations (21). 

If, now, a further term is added to O, the equations for the 


variations of the arbitrary constants are (by [29]) seen to be 



which form a canonical set. 

In particular, if C„ +1 , C„ + 2 , ■ . • C 2)t , are all zero or definite 
constants, we see, by putting 



in which 4>=o, that a„ a 2 , . . . a,., C 2 , . . . C„, in this case form 

a canonical set of variables. For example, in the case of the 
planetary theory, 


a, 9, 


a 



form a canonical set of arbitrary constants of the elliptic motion. 
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Note on Mr. A. T. G. Campbell’s paper “ On the Variation of 
TJncanonical Arbitrary Constants, with an Application to the 
Planetary Theory.” By Sir Robert S. Ball, LL.D., P.R.S. 


The object of this note is to set down another method by 
which the calculation of Lagrange’s “parentheses” in the 
planetary theory can be effected very simply. The method was 
suggested by a manuscript which Mr. Campbell kindly com¬ 
municated to me. I have also to acknowledge the assistance 
derived from the instructive remarks of Mr. Cowell and Mr. 
Whittaker, to whom, as well as to Mr. Campbell, it has been my 
privilege to have lectured at Cambridge. 

Adopting the notation of Mr. Campbell’s manuscript, we put 

0 (x, £•) ^ dx dx dx dx 
0 (A., fx) dx du. du. d\’ 


where X, p are any two of the six arbitrary constants. 

We thus have the expression for the “ parenthesis ” 

rx in - i ) + v) + i ) 

L ’ J 

We now make the transformation to another set of rect¬ 
angular coordinates x x , y x , z x , where 

x — l\X, 

y = m l x 1 + m 2 y l + m 3 ~ l 

s — + n.,y, + >i 3 c l 

It is unnecessary to write down the general result. For the 
planetary theory we take the special case where 


l j= cos a 

l 2 = — sin a cos <p 

£3 — 0 

m x = sin a 

m 2 — cos a cos <p 

™ 3 =° 

n x — 0 

= sin <p 

?i z = O 


And if we express x x by £ and y x by y we get after a little 
reduction 


[A, #0 = 


0(1, $) + d(iM» 

0 (A) M) 0 (A, M) 


- sin <p{£fi-in) 


0(«, ft) 

0(A, /*) 


da d / j. • ji \ . da d s >. > \ 

+ cos * S. T,. W-“ s f ^ » 


We now transform to two new variables X and Y given by 
the equations 

£ = X cos g — Y sin g 
rj = X sin g -f Y cos g 
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